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Abstract 

. High-precision measurements of the relic dark matter density and the calculation of dark matter 

> . 

f — , annihilation branching fractions in the sun or the galactic halo today motivate the computation of 

CN ' the neutralino annihilation cross section beyond leading order. We consider neutralino annihilation 



via squark exchange and parameterize the effective annihilation vertex as a dimension-six operator 
suppressed by two powers of the squark mass and related to the divergence of the axial vector 
Q_i| current of the final-state quarks. Since the axial vector current is conserved at tree level in the 



limit of massless quarks, this dimension-six operator contains a suppression by the quark mass. The 
quark mass suppression can be lifted in two ways: (1) by corrections to the dimension-six operator 



involving the anomalous triangle diagram, and (2) by going to dimension-eight. We address the first 
of these possibilities by evaluating the anomalous triangle diagram, which contributes to neutralino 
annihilation to gluon pairs. We relate the triangle diagram via the anomaly equation to the decay 
of a pseudoscalar into two gluons and use the Adler-Bardeen theorem to extract the next-to-leading 
order (NLO) QCD corrections to XX ~^ 99 from the known corrections to pseudoscalar decay. The 
strong dependence of the dominant XX ~^ QQ cross section on the relative velocity of the neutralinos 
makes these NLO corrections unimportant at x decoupling but significant today. 
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I. INTRODUCTION 



The presence of non-baryonic dark matter in the universe is compelling evidence for 
physics beyond the Standard Model. While many very different models have been proposed 
to explain the dark matter (wimps, axions, wimpzillas, modified gravity, etc.) the 
thermal production of stable weakly-interacting particles with weak-scale mass remains as 
the most attractive and predictive explanation for the observed dark matter relic abundance, 
and further allows the solution to the dark matter problem to be linked to the solution to the 
hierarchy problem and tested at current and future collider experiments. Supersymmetry 
provides an especially attractive explanation with the lightest super symmetric particle (LSP) 
as the dark matter candidate. Once such a weakly-interacting massive particle (x) has 
been discovered and its couplings measured, it will be possible to compute its annihilation 
cross section (which controls its thermal relic abundance) and compare to the measured 
dark matter abundance to test our understanding of the microphysics of dark matter. The 
cosmological dark matter abundance is already measured at the 10% level to be Q] 

fi CDM /i 2 = 0.12 ± 0.01 (SDSS+WMAP) 

and future cosmic microwave background experiments such as PLANCK expect to improve 
this to the few-percent level |2j. In order to match the expected few-percent precision of the 
cosmological measurements, we need both high-precision inputs from the colliders and high- 
precision calculations of the neutralino annihilation cross section. The latter requirement 
means going beyond leading order. High-precision calculations of the relic abundance are 
thus needed to match the microphysics onto cosmology. 

Another important physics application of these higher order QCD corrections is the cal- 
culation of signals from WIMP annihilations in the galatic halo or in the interior of the sun. 
The NLO corrections are potentially important in the evaluation of the branching fractions 
for the observable gamma ray and neutrino signals, respectively. As we shall see, these 
NLO corrections turn out to be much more important for the calculation of the observable 
gamma ray and neutrino signals than in the calculation of the relic density because of the 
strong dependence of the tree-level annihilation cross-section on the relative velocity of the 
neutralinos. 
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FIG. 1: Feynman diagrams that contribute to the total annihilation cross section: (a) the tree level 
diagram, (b-d) t-channel squark exchange, and (e-j) s-channel Z and Higgs exchanges. 

A. Neutralino annihilation cross section 

The behavior of the annihilation cross section depends on the composition of the neu- 
tralino. Throughout this paper we assume that the LSP is largely gaugino as motivated by 
mSUGRA models jj| . The processes that contribute to the cross section up to order o? s and 
one loop are shown in Fig. The tree- level diagram is shown in Fig. Q^a). Figures ^b-d) 
show the diagrams with t-channel squark exchange, whereas (e-j) show the diagrams with 
s-channel Z, H°, h°, A exchanges. The gauge and Higgs bosons couple to the Higgsino part 
of the LSP and thus their contributions are suppressed for a mostly-gaugino neutralino. 1 
The corresponding suppression factors for the s- and p-wave terms in the cross section are 
given in Table UJ 

B. The anomaly equation 

The leading contribution to neutralino annihilation via exchange of a squark of mass 
M, shown in Fig. HJa), can be reduced to an effective vertex described by a dimension-six 
operator suppressed by M 2 , 

C = (c/M 2 )0 6 , O e = (xi»lsX)(qYlsq) , (1) 
1 The Higgsino fraction suppression can be removed at the cost of going to one loop jfj. 
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TABLE I: Dependence of the cross section from each diagram on various suppression factors. The 
diagrams are shown in Fig. ^ The p-wave contribution to the cross section is shown only when 
it is important due to suppression or absence of the s-wave component. The column labeled "x 
mixing" sketches the dependence of the cross section on the neutralino composition. Interference 
terms between the various diagrams carry a combination of the suppression factors corresponding 
to each diagram, and are not shown. We note that the Z-pole structure of diagram (i) is canceled 
by a numerator factor supplied in accordance with Yang's theorem P, Q]. 

where c is a dimensionless coefficient. This dimension-six operator corresponds to taking the 
leading term in the expansion of the squark propagator in powers of I/M 2 ; in particular, 
we work in the limit m 2 -C M 2 . 

In the static limit, where the relative velocity of the two neutralinos can be neglected, 
the operator Oq is related to the divergence of the axial vector current of the quarks qq: 



(2) 



In the massless quark limit, m q = 0, the axial vector current is conserved at tree level, 
^(97^759) = 0, and all tree amplitudes due to the dimension-six operator vanish; in par- 
ticular, radiating additional gluons cannot lift the suppression. Even at the loop level, for 
example diagrams involving the exchange of a virtual gluon, the suppression is still valid un- 
less the anomalous triangle diagram is involved. This is the well-known partially-conserved 
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axial current (PCAC) condition. 

Indeed, only through the anomalous loop diagrams is the conservation of the axial vector 
current violated, in the form ^ 

d,(qYl 5 q) = 2m q qz l5 q + ^-G^&^ (3) 

with \G iiV = e^ va pG a P denoting the dual color field strength tensor. The simplest such 
anomalous diagram is shown in Fig. H^c). Neglecting the mass of the internal quark q' and 
using the anomaly equation, this diagram can be written in the form 

Ce,(xX - 99) = (^) (X *75 X)|GWGW^ , (4) 

for m q i <C m x . In the opposite limit, m q i ^> m x , the very heavy quark decouples; the top 
quark contribution can be neglected if m x < 100 GeV. The leading order (LO) calculation 
using the anomaly equation was first studied correctly by Ref. |Q| in the 77 channel in QED. 

The gluonic decay amplitude of a fundamental pseudoscalar, A — > gg, is also related 
to the anomaly equation. This decay proceeds through a quark loop. In the heavy quark 
limit, ttlq ^> m^, the divergence term on the left-hand side of Eq. Q becomes insignificant, 
leading to, 

~ 2m Q Qz l5 Q + SjLGMGMr (5) 

Note that because the Yukawa coupling of A to the quark Q is proportional to the quark 
mass, the tuq dependence of the A — > gg partial width drops out in the limit vhq 3> vn^. 

In contrast, the neutralino pair annihilates into two gluons via the anomaly diagram, 
which is dominated by light quarks, m q -C m x . In this limit, the term proportional to m q 
on the right-hand side of Eq. becomes insignificant, leading to, 

d.qYlM * + ^-G^G^. (6) 

Thus we see that two seemingly different processes, A — > gg through a heavy quark loop 
and XX ~ *■ 99 through a light quark loop, are related by the anomaly equation to the same 
gluonic operator. The Adler-Bardeen theorem js| guarantees that the anomaly equation, 
Eq. (jHJ), is valid to all orders of a s . One can take advantage of this anomaly property to 
obtain the higher-order QCD corrections to XX ~ > 99 from the known results for A — ► gg 
at next-to-leading order (NLO) from Ref. |10] . Note that, because of the non-abelian nature 
of the gauge field, the above gluonic operator also incorporates tri-gluon amplitudes beyond 
leading order. 



C. Beyond dimension-six 



If we include higher terms in the 1/M 2 expansion, the PCAC constraint will be lifted. 
The dimension-eight operator corresponding to an amplitude proportional to 1/M 4 survives 
even in the massless quark limit, m q — > 0. One can use this dimension-eight amplitude 
to calculate the rate of XX ~ * QQ9 from diagrams such as Fig. ^b); a full calculation was 



done in Ref. 11 1. However, the contribution to XX QQ9 from the dimension-six operator 
due to the anomaly with a virtual gluon turning into a quark pair [Fig. W[c)} suffers less M 
suppression and will dominate the dimension-eight term for M 3> m x even though the order 
in a s is higher. 

The effective vertex of the dimension-eight operator for XX QQ9 is 

£ = (c 8 g 2 s /M 4 )0 8 , 8 = e, ua(3 G^(qYl 5 T c q)( X Yl 5 x), (7) 

where G^ al3 is the gluonic field strength for color index c and T c is the corresponding SU(3) 
generator. This operator has exactly the same form as that in Eq. (6) of Ref. [l2j for the 
XX ~^ ffl amplitude computed through explicit expansion of the propagators to order 
1/M 4 in the limit m q = 0. 

It is interesting to note that the amplitude for XX ~^ QQ9 from the dimension-six operator 
with one gluon splitting into qq [shown in Fig. ^d)] yields an operator of the same form as 
in Eq. (0). This allows the interference term between this diagram and the dimension-eight 
process to be easily obtained. 

II. CALCULATION 

A. The anomaly at leading order 

Since neutralinos are Majorana in nature, the initial state behaves as a pseudoscalar in 
the zero- velocity limit Q| . In particular, for v Ie \ = the antisymmetrized neutralino spinors 



reduce to the projection operators |l z 

u(pi)v(p 2 ) - u(p 2 )v(p 1 ) = (m x + f )7 5 = m x (l + 7°)7 5 . (8) 
where p\ and p 2 are the four- momenta of the incoming neutralinos and 2P = p\ + p 2 . 
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We work in the limit of zero fermion mass, in which case the off-diagonal terms in the 
squark mass matrices vanish and the squark mass eigenstates coincide with the electroweak 
eigenstates q^ and (Jr. Applying the reduction formula Eq. JSj) to the amplitude of the 
XX 99 diagram shown in Fig. ^d) allows us to write the amplitude for the diagram 
involving squark <^ as 

-^9 2 s9 2 r ,i f d 4 q 



Ml 



J Wji Tl [^)r ai T75^v] <h);e(k 2 )U (9) 



where the neutralino-quark-squark couplings are defined for right- and left-handed squarks 
respectively as 



9r 



V2AWQ, 91 = -V2N n g'{T 3 - Q) + V2T 3 N 12 g. (10) 



Here T 3 is the squark isospin, Q is the squark electric charge, and N u and N 12 are the 
bino and wino components of the neutralino as defined in Ref. [l5|. We also define Pr/l — 
(1 ±7 5 )/2 in the usual way as the right- and left-handed projection operators. The external 
gluon momenta are called k\ and k 2 , and q is the momentum flowing in the loop. 

The form factor J r ( y q)^ u,ab for squark qi and the corresponding internal quark g« is given 
explicitly by 



<j- ¥2 + m gt b q 1 + m qt ^ a 4 + ¥1 + m gi 
(q — k 2 ) 2 — ml. q 2 — ml. (q + ki) 2 — m 2 



(/-¥i + m qi r u f g ^+ m qi b g'+¥2 + m qi 



(q — k\) 2 — m 2 . q 2 — m 2 (q + k 2 ) 2 — m 2 

The two terms in J-(q)f u ' a correspond to the two possible directions of fermion flow in 
Fig. IHd). The diagrams with the neutralinos crossed are already included through the use 
of the antisymmetrized spinors in Eq. (jSJ). 

After summing over left and right squark states, the amplitude becomes 

Mt = -J ^2Tr [f(C'V(V, + An,)] PAh)l<k 2 )l, (12) 

with 

v _V2gl(jl gf\ , V2 9 2 (g 2 gf\ 

Vl ~ 4 \M 2 M 2 ) Ai ~ 4 \M 2 + M 2 ' 

\ Qi r lil J \ qir q%l / 

with couplings g r j given in Eq. ([TU]l. 

After integrating Eq. (JT^J), the piece involving the vector coupling vanishes due to the 
conservation of vectorial current (CVC): the contracted vector behaves as a divergence, so 



that the resulting vector coupling to multi-gluon states vanishes. Light quark masses can be 
neglected in the form factor if ml -C m 2 and the top quark loop amplitude is suppressed if 
vr? x <C mf . In these limits, the loop amplitude sums over five massless quarks. Note that in 
the massless quark limit, the form factor J r (q)^' u,ab becomes independent of the quark flavor 
i. The sum over quarks in the loop can then be factorized into a sum over the coupling 
factors Ai times the universal massless form factor J r (q)f u,ab . 

The gluon production amplitude contains a pseudovector triangle diagram. This can be 
transformed to a pseudoscalar triangle diagram via the axial anomaly, Eq. (J3J), where = 
—2iP^. The anomaly is computed by relating it to the decay of a fundamental pseudoscalar 
Higgs boson A to two gluons via a heavy quark loop. If the mass of the heavy quark Q is 
sufficiently large, tuq 3> m^o, then from Eqs. (0) and (jUJ) we have, 

2iP a A iQa a l5Qi = 2m Q Qi^ 5 Q ^ A i- ( 14 ) 

The amplitude becomes 

^Co(0, 0, s, m% m 2 Q , m 2 Q )Tr [t a t b ] e^k^k^k^k^ (15) 

where s = 4m£ and Cq(0, 0, s, mg, mg, m^) is a three-point Passarino-Veltman integral (if] ]. 
In the limit of heavy quark mass, 3> s, the three-point integral reduces to —l/2m 2 Q. In 
this limit the dependence on the heavy quark mass drops out and the amplitude becomes 

M = V^e^% a ^e(fc0;e(fc 2 ):, (16) 

where we have used Tr[t a t b ] = (l/2)5 ab . Squaring the amplitude and integrating over phase 
space gives the leading-order (xX ~^ 99) annihilation cross section, 

jjl^ / \ 

Vrcicr lo(xX -> 99) = \^ Ai J ' (17) 

where in our approximation the sum runs over the five light quarks gj. Our result agrees 
with that of, e.g., Ref. in the limit m q = 0. 
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B. Beyond leading order 

As discussed in Sec. II Bl we can use the Adler-Bardeen theorem Jj| and Eqs. (0) and © 
to obtain the QCD corrections to the XX 99 annihilation cross section by exploiting the 
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known results for pseudoscalar Higgs decays to gluon pairs, A — > gg, beyond leading order. 

The NLO QCD corrections to the A — > gg partial width were calculated by Spira et 
al. 1(3 ]. In the heavy top quark limit, for which our anomaly relation is valid, the NLO 
QCD corrections are given by a multiplicative factor |lOj ]. times the LO decay rate 



T N lo(A° -> gg) = T LO (A° 



1 + 



7T 



97 7 Ar 33-2iV fl u 2 
4 6 6 m \ j j 



(18) 



where /x is the renormalization scale. The integer Nf counts the number of quark flavors 
in the gluon splitting, with Nf = 5 for m& m x m t . The diagrams that contribute to 
A —* gg at LO and NLO are shown in Fig. |21 The NLO final states include gg, ggg, and 

gqq- 



A° 





ifo7Rr> g -. 






(f) 




< 



(h) 







FIG. 2: LO and NLO diagrams for A — ► 55. (a) is the leading order process, (b-d) are the real 
emission diagrams with three final-state particles, and (e-1) are virtual corrections to diagram (a). 



In XX ~ > 99 & t NLO, a divergence occurs for the diagram in Fig.^d) when the final-state 
quarks are soft or collinear, in which case the gluon propagator diverges. This is the source 
of the logarithmic enhancement factor, log(m?./m^), found for this diagram in Ref. 



121. 



However, this logarithmic term is precisely canceled by the renormalization of the strong 




FIG. 3: XX ~> 99 diagrams that supply a logarithmic factor to cancel that from Fig. ^ The 
one-loop corrections to the gluon legs include quark and gluon loops. 

coupling due to the quark bubble that appears in the virtual part of the NLO correction, 
shown as the interference of the diagrams in Fig. El This is the familiar cancellation of 
logarithmic divergences guaranteed by the Kinoshita-Lee-Nauenberg theorem Q). A similar 
cancellation occurs for the analogous diagrams in which the soft /collinear quarks are replaced 
with gluons. 

We now invoke the Adler-Bardeen theorem [9!] and take over the NLO corrections to 
^4° ~~ * 99 to the XX ~ > 99 process in the zero- velocity limit. In this correspondence, the 
pseudoscalar mass tua is replaced by the XX center-of-mass energy, equal to 2m x in the zero- 
velocity limit. The NLO correction to the cross section for XX ~^ 99 follows immediately 
from Eqs. (fT7jl and (fT%j) . 

m l „\ 2 L « s /97 7 Ar 33-2iV fl /x 2 \1 . . 

We give explicitly the result for fi = 2m x and Nf = 5, 

VrelCTNLoiXX ^99) = ^ (1 + 0.62), (20) 

where we set m x = 100 GeV. The strong coupling is evaluated based on five-flavor running 
at the scale \x = 2m x where it appears both explicitly and within the coefficient Ai. We note 
that the above choice of m x = 100 GeV is well above the current experimental limit [rsj ]. 
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III. NUMERICAL RESULTS 



In this section we examine the validity of our assumption of massless quarks in the loop for 
the LO XX ~ > 99 calculation, show the improvement in the renormalization scale dependence 
obtained in going to NLO, and compare the XX ~^ 99 annihilation cross section to that of 
the leading-order tree level process XX ~^ through t-channel squark exchange. 

We have assumed in our use of the anomaly equation that the five light quarks running 
in the loop for XX ~^ 99 were massless, and we neglected the heavy top quark contribution. 
In Fig. El we test this assumption for the LO cross section by comparing our approximation 
to the full cross section including the quark mass dependence (we continue to use the 1/M 2 
approximation for the squark propagator). We plot the full cross section normalized to our 
five-massless- quark approximation as a function of m x . The full formula differs from our 
approximation by less than 10% for 6 GeV < m x < 110 GeV. For heavier neutralinos, the 
top quark loop starts to have a significant effect, with destructive interference occurring 
between the top loop and the lighter quark loops. Two of the top quarks in the loop go 
on shell at m x = m t , leading to the large dip in the cross section. For neutralinos lighter 
than about 6 GeV, the nonzero mass of the bottom quark begins to play a significant role, 
leading to the dip at lower masses. 




50 100 150 200 



m z (GeV) 

FIG. 4: The effect of quark masses on v re iaLo(xX ~* 99)- Shown is the full cross section including 
the quark mass dependence, normalized to the approximate cross section obtained by setting 
m q = for the five light quarks and mt = oo, for a common squark mass M = 200 GeV and a 
pure bino neutralino, Nu = 1, N%j = (J ^ 1). 
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In an all-orders calculation, physical observables cannot depend on the renormalization 
scale fj,. The \i dependence of our predictions is an artifact of computing to a finite order 
in perturbation theory. The \x dependence can then be used to estimate the size of the 
uncomputed higher-order corrections. The dependence of the XX ~ * 99 annihilation cross 
section on the renormalization scale at LO and NLO is shown in Fig. E(a) for M = 200 
GeV, m x = 100 GeV, and a pure bino neutralino, Nn = 1, Nij = (j ^ l). 2 

Varying \x by a factor of two in either direction from the central value \x = %rn x yields a 
scale dependence of ±16% at LO and ±9% at NLO. The corresponding annihilation cross 
sections including scale uncertainty are shown in Fig. Efb) as a function of m x . 
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FIG. 5: Cross section for XX ~> 99 a t LO and NLO for a pure bino neutralino with M = 200 GeV. 
(a) Dependence on the renormalization scale \i for m x = 100 GeV. (b) Dependence on m x showing 
the renormalization scale dependence in the band \x = m x , 4m x . 



Exact cross section formulae for all tree-level two-to-two neutralino annihilation processes 
are given in Ref. Q|. Using the expansion of the thermally averaged cross section in terms 
of x — T/m x , 

{v re ia) = a + bx, (21) 

one can compare the leading tree-level process, XX ~^ QQ y i a ^-channel squark exchange, to 
our results for XX ~^ 99 &1 NLO given in Eq. (j2*U|) both for annihilation during freeze-out 

2 We note here that the renormalization scale dependence of the tree-level XX ~ * cross section arises only 
from the running quark mass in the s-wave contribution Since the cross section for this process is 
dominated by the quark-mass-independent p-wave part during freeze-out in the early universe, the scale 
dependence is negligible at tree-level. 
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in the early universe, x ~ 1/20, and for annihilation in the galactic halo today, x ~ 
(corresponding to v/c ~ 10~ 3 ). The cross sections for the two annihilation processes are 
shown in Fig. EJa) as a function of m x for M = 200 GeV and a pure bino neutralino, with 
x = 1/20 and x = 0. For XX ~ y QQ we sum t ne cross section over the five light final-state 
quark flavors and neglect XX annihilation into lepton pairs. We use the running quark mass 
m q (fi), which serves to resum the leading logarithmic QCD corrections to this process from 
soft gluon radiation [llj], and take the renormalization scale fi = 2m x . The annihilation 
cross section for XX ~^ 99 is dominated by the s-channel component, and so we show only 
one curve for x — 1/20 and x — 0. The tree- level cross section for XX ~^ QQ depends strongly 
on the relative velocity of the neutralinos, because the s-channel cross section is suppressed 
by the final-state quark mass. Thus the annihilation cross section at x — 0, which comes 
only from the s-wave component, is quite small and is comparable to that from XX ~~ *■ 99- 
At x = 1/20, on the other hand, the XX ~^ cross section is dominated by the p-wave 
component and is larger than XX ~^ 99 by almost a factor of 100. 

In Fig. Efb) we show the corresponding K-factors, defined as the ratio of the total cross 
section, XX ~ > <?<? + 99-, to the tree-level XX ~^ <?<? cross section. This gives a measure of 
the relative importance of the XX 99 component of the total annihilation cross section. 
We see that during freeze-out, x ~ 1/20, the XX 99 contribution is quite small and K = 
1.01 — 1.02 over the range of m x considered. In the present epoch, however, with x ~ 0, the 
K-factor is considerably larger, K = 1.3 — 30 depending on the mass of the neutralino. Such 
a large K-factor will impact annihilation branching fractions today, changingthe gamma 
ray flux from the galactic halo and the neutrino flux from inside the Sun |l|, lll| . We note 
also that because both the XX ~^ QQ an d XX ~^ 99 annihilation cross sections have the same 
leading 1/M 4 dependence on the squark mass, these K-factors will not depend significantly 
on the common squark mass scale. 

We now exhibit the dependence of the annihilation cross section on the neutralino compo- 
sition. Because we have worked in the zero-quark-mass limit in our calculation of XX ~^ 99 S 
we have neglected the couplings of Higgsinos to the internal quark loop, which are propor- 
tional to the quark mass. We thus consider only mixed bino-wino neutralinos. We can then 
parameterize the mixing coefficients Ny in terms of a bino-wino mixing angle as 

iVn = cos 0, iVis = sin 0, iV 13 = iV 14 = 0. (22) 
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25 50 75 100 ' 20 40 60 80 100 



m % (GeV) m z (GeV) 

FIG. 6: (a) Annihilation cross sections for XX ~^ QQ through i-channel squark exchange [dia- 
gram IB^a)] and XX ~* 99 a t NLO for M = 200 GeV and a pure bino neutralino. We show XX ~^ QQ 
for both x = 1/20, corresponding to freeze-out in the early universe, and x = 0, corresponding to 
annihilation in the galactic halo at the present time. The steps in the XX ~~ QQ cross section for 
x = at low m x are due to quark mass thresholds, (b) Corresponding K- factors. We plot K — 1, 
which is the ratio of annihilation cross sections of XX ~^ 99 to XX ~^ QQi f° r x = and x = 1/20 as 
shown in (a). 



In Fig. UJ we again compare the leading tree-level process, XX ~^ QQ v i a ^-channel squark 
exchange, to our results for XX 99 a t NLO. The cross sections for the two annihilation 
processes are shown in Fig. Ufa) as a function of the bino-wino mixing angle </> for m x = 
50 GeV, and M = 200 GeV, with x = 1/20 and x — 0. Again we take the renormalization 
scale fi = 2m x . There is a large enhancement of both annihilation cross sections if the 
neutralino has a large wino component, ~ 90°, due to the stronger coupling of the wino 
to a quark-squark pair. In Fig. [7(b) we show the corresponding K-factors. The nontrivial 
dependence of the K-factors on arises from the interference among the five light quark 
loop diagrams that contribute to XX ~^ 99- By contrast, there is no interference between the 
amplitudes for XX ~ y QQ with different flavor squarks in the t-channel because the internal 
squark flavor is fixed by the external quark flavor. 

Finally, we note that the next-to-next-to-leading order (NNLO) corrections to A — ► gg 
have been computed in Ref. 20(. One may be tempted to take this correction over to XX ~^ 
gg in the same way as the NLO correction. However, at NNLO the A — > gg decay receives a 
contribution from the interference between the G^}G^ a ^ u operator and a d^qj^^q operator 
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FIG. 7: (a) Annihilation cross sections for XX ~^ QQ an d XX ~> 99 as m Fig. as a function of the 
bino-wino mixing angle (j), for m x = 50 GeV. Pure bino corresponds to <p = 0°, 180° and pure wino 
to 4> = 90°. (b) Corresponding K- factors. 



generated at two-loop level. Once effective operators other than G ( j^G^ fJ-l/ appear, our use 
of the anomaly equation no longer applies. A proper treatment of XX ~ > 99 at NNLO would 
thus require a new calculation of the operator matching conditions and renormalization. 



IV. CONCLUSIONS 



We reviewed the dependence of the main neutralino annihilation processes on various 
suppression factors - the Higgsino fraction, the quark and squark masses, and the relative 
neutralino velocity - and identified the dominant XX ~~ QQ annihilation process for a gaugino- 
like neutralino as due to a dimension-six operator in the zero- velocity limit. This dimension- 
six operator contains the divergence of the axial vector current of the quarks qq, which leads 
to the well-known quark mass suppression of the annihilation cross section. This quark 
mass suppression can be lifted in two ways: either through corrections to the dimension- 
six operator involving the anomalous triangle diagram, or by going to dimension-eight. We 
focused on the anomalous triangle diagram, which describes neutralino annihilation to gluon 
pairs. In the approximation of massless quarks running in the loop, the anomaly equation 
relates XX ~ > 99 t° the seemingly unrelated process of pseudoscalar decay to gluon pairs 
via a very heavy quark loop. We used this relation to compute XX ~^ 99 m terms of the 
decay process A — > gg. Further, taking advantage of the Adler-Bardeen theorem which 
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guarantees that the anomaly equation is valid to all orders in a s , we extracted the NLO 
QCD corrections to XX 99 from the known corresponding results for A — ► gg and 
wrote them as a simple multiplicative factor that can be easily inserted into numerical 
neutralino annihilation codes. For m x = 100 GeV and fi = 2m x , the NLO QCD corrections 
increase the annihilation cross section by 62%. The NLO corrections also reduce the residual 
renormalization scale dependence of the XX ~^ 99 annihilation cross section from ±16% to 
±9%. 

Our NLO results were computed in the approximation <C m x <C m t . This ap- 
proximation yields a LO XX ~~ * 99 cross section within 10% of the exact result for 
6 GeV < m x < 110 GeV. We finally compared our results for XX ~^ 99 a t NLO to the 
dominant XX QQ cross section in this neutralino mass range. For neutralino annihilation 
during freeze-out in the early universe, our results for XX ~^ 99 & t NLO constitute only 
1 — 2% of the dominant cross section for a gaugino-like neutralino and are thus of little 
importance for computing the relic neutralino abundance. However, for neutralino annihi- 
lation at the present time the relative neutralino velocity is much lower, leading to a much 
smaller tree-level XX ~^ <?<? cross section. In this situation, the XX ~^ 99 cross section can 
be as large or larger than XX ~^ QQi so that NLO corrections can have a significant impact 
on the computation of gamma ray and neutrino fluxes from neutralino annihilation in the 
galatic halo and inside the Sun, respectively. 
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